OU-TAP 115: ApJ in press 



The Variation of Gas Mass Distribution in Galaxy Clusters: Effects of 

Preheating and Shocks 

Yutaka Fujita 1 and Fumio Takahara 

Department of Earth and Space Science, Graduate School of Science, Osaka University, 
q \ Machikaneyama-cho, Toyonaka, Osaka, 560-0043 

(N ; 

a5 1 ABSTRACT 

We investigate the origin of the variation of the gas mass fraction in the core of 
galaxy clusters, which was indicated by our work on the X-ray fundamental plane. 

■ Applying a spherical collapse model of cluster formation and considering the effect of 

shocks on preheated intracluster gas, we construct a simple model to predict the spatial 
^ . gas distribution of clusters. As is suggested by our previous work, we assume that the 

core structure of clusters determined at the cluster collapse has not been much changed 
after that. The adopted model supposes that the gas distribution characterized by the 

\ slope parameter is related to the preheated temperature. Comparison with observations 

Q-f of relatively hot (> 3 keV) and low redshift clusters suggests that the preheated temper- 

ature is about 0.5-2 keV, which is higher than expected from the conventional galactic 
wind model and possibly suggests the need for additional heating such as quasars or 
gravitational heating on the largest scales at high redshift. The dispersion of the pre- 

. heated temperature may be attributed to the gravitational heating in subclusters. We 

calculate the central gas fraction of a cluster from the gas distribution, assuming that 
the global gas mass fraction is constant within a virial radius at the time of the cluster 
collapse. We find that the central gas density thus calculated is in good agreement with 
the observed one, which suggests that the variation of gas mass fraction in cluster cores 
appears to be explained by breaking the self-similarity in clusters due to preheated gas. 
We also find that this model does not change major conclusions on the fundamental 
plane and its cosmological implications obtained in previous papers, which strongly sug- 
gests that not only for the dark halo but also for the intracluster gas the core structure 
preserves information about the cluster formation. 
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1. Introduction 



Correlations among physical quantities of clusters of galaxies are very useful tools for studying 
the formation of clusters and cosmological parameters. Recently, we have found that clusters at low 
redshifts (z < 0.1) form a plane (the fundamental plane) in the three dimensional space represented 
by their core structures, that is, the central gas density /9 gas ,cb core radius r c , and X-ray temperature 
Tgas (Fujita and Takahara 1999a, Paper I). On the other hand, a simple theoretical model of cluster 
formation predicts that clusters should be characterized by the virial density p v [ r (or the collapse 
redshift z co ii) and the virial mass M W - 1T (Fujita and Takahara 1999b, Paper II). Thus, assuming the 
similarity of the dark matter distributions, clusters should form a plane in the three dimensional 
space of the dark matter density in the core pdm.c, the core radius of dark matter distribution 
rDM,d and the virial temperature T v ; r 2 . However, the relations between the two planes are not 
simple; for example, it is found that p gas ,o is not proportional to Pdm,c- in Paper I, we found that 
the ratio p gas ,o/PDM,c is not constant but obeys the relation of /O ga s,o/PDM,c oc PdmV^dmc' where 
Mdm, c is the core mass. This raises the question how the segregation between gas and dark matter 
occurs. 

In the hierarchical structure formation, dark halos are expected to obey scaling relations. In 
fact, numerical simulations suggest that the density distribution in dark halos take a universal 
form as claimed by Navarro et al. (1996, 1997). On a cluster scale, it can be approximated by 
PdmM of r ~ 2 for r < 1 Mpc, where detailed X-ray observations have been done (Makino et al. 
1998). On the contrary, observations show that the slope of the density profile of the hot diffuse 
intracluster gas has a range of value. Radial surface brightness profiles of X-ray emission are often 
fitted with the conventional (5 model as 

I{R) = (1+^/^)3^-1/2 ' (1) 

where /3 Q b s is the slope parameter (Cavaliere and Fusco-Femiano 1978). If the intracluster gas is 
isothermal, equation (1) corresponds to the gas density profile of 



Pgas(r) ~ ( 1+r 2/ r 2)3/W2 • (2) 

Observations show that the slope parameter takes a range /3 b s ~ 0.4 — 1 (Jones and Forman 1984, 
1999). This means that for r >> r c , the density profiles range from oc r~ 1,2 to oc r~ 3 , which are 
more diverse than those of dark matter. Moreover, observations show that the clusters with large 
r c and T gas tend to have large /3 Q bs (e.g. Neumann and Arnaud 1999; Horner et al. 1999; Jones and 



2 In Paper I, we used the terms ' virial density', 'virial radius', and ' virial mass' to denote the dark matter density 
in the core, the core radius of dark matter distribution, and the core mass, respectively. This is because we assumed 
that the dark matter density in the core is proportional to the average dark matter density over the whole cluster 
(Paper II). To avoid possible confusions, here we use the term ' dark matter', and the term ' virial' will be used to 
represent spatially averaged quantities of gravitational matter (mostly dark matter) within the virialized region. 
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Forman 1999). Since the average gas fraction of clusters within radii much larger than r c should 
be universal and the dark matter distribution of clusters is also universal, the variation of /? t, s 
is expected to correlate with that of the gas fraction in the core region. In other words, the gas 
fraction at the core is not the same as that of the whole cluster and is not proportional to the dark 
matter density. This fact must be taken care of when we discuss cosmological parameters using 
observational X-ray data. Since the emissivity of X-ray gas is proportional to p gas , most of the 
X-ray emission of a cluster comes form the central region where p gas is large. Although in Papers 
I and II, we did not take account of the effects of /3 bs) we did find the gas mass fraction in the 
core region is diverse by analyzing the X-ray emission. In this paper, we reanalyze the data taking 
account of /3 Q b s and discuss the relation between core and global gas mass fractions. We will also 
show that major conclusions on the fundamental relations are not changed. 

The variation of gas mass fraction itself has been investigated by several authors (e.g. Ettori 
and Fabian 1999; Arnaud and Evrard 1999). Ettori and Fabian (1999) argue that it is partially 
explained if the dark matter has a significant baryonic component. Another possible explanation of 
the diverse gas distributions is that intracluster gas had already been heated before the collapse into 
the cluster; the energetic winds generated by supernovae are one possible mechanism to increase gas 
entropy (e.g. Dekel and Silk 1986; Mihara and Takahara 1994). In fact, Ponman et al. (1999) find 
that the entropy of the intracluster gas near the center of clusters is higher than can be explained 
by gravitational collapse alone. In order to estimate the effect of the preheating on intracluster 
gas, we must take account of shocks forming when the gas collapses into the cluster; they supply 
additional entropy to the gas. Cavaliere et al. (1997, 1998, 1999) have investigated both the effects 
and predicted the relation between X-ray luminosities and temperatures (Lx — T gas relation) . They 
predicted that the gas distributions of poor clusters are flatter than those of rich clusters, which 
results in a steeper slope of Lx — T gas relation for poor clusters. This is generally consistent with 
the observations. It is an interesting issue to investigate whether this scenario provides a natural 
explanation for the observed dispersion of gas mass fraction in the cluster core and whether it 
reproduces the X-ray fundamental plane we found in Paper I in our general theoretical scenario. 

In order to clarify what determines the gas distribution, we construct as a simple model as 
possible. Although many authors have studied the preheating of clusters (Kaiser 1991; Evrard 
and Henry 1991; Metzler and Evrard 1994; Balogh et al. 1999; Kay and Bower 1999; Wu et al. 
1999; Valageas and Silk 1999), this is the first time to consider the influence of the preheating and 
shocks on the fundamental plane and two-parameter family nature of clusters paying attention to 
the difference between the collapse redshift z co \\ and the observed redshift z \^ s of clusters explicitly. 
In §2, we explain the model of dark matter potential and shock heating of intracluster gas. In §3, 
we use the model to predict (3 Q ^ S — T gas and (3 Q ^ S — r c relations, and the fundamental plane and 
band of clusters. The predictions are compared with observations. 
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2. Models 

2.1. Dark Matter Potential 

In order to predict the relations among parameters describing a dark matter potential, we use 
a spherical collapse model (Tomita 1969; Gunn and Gott 1972). Although the details are described 
in Paper II, we summarize them here for convenience. 

The virial density of a cluster p YU at the time of the cluster collapse {z co \\) is A c times the 
critical density of a universe at z = z co \\. It is given by 

A I \ A Z7V \ 2 a ^o(l + ^coll) 3 / q \ 

Pvir — ^cPcritl^collJ — ^ePcrit,0- fc \ z co\\ ) ~ ^ePcrit.O ^ , l<jj 

where Q(z) is the cosmological density parameter, and E(z) 2 = f2o(l + z) 3 /Q(z), where we do not 
take account of the cosmological constant. The index refers to the values at z = 0. Note that the 
redshift-dependent Hubble constant can be written as H(z) = W0hE(z) km s" 1 Mpc -1 . We adopt 
h = 0.5 for numerical values. In practice, we use the fitting formula of Bryan and Norman (1998) 
for the virial density: 

A c = 18vr 2 + 60x - 32x 2 , (4) 

where x = fi(z co n) — 1. 

It is convenient to relate the collapse time in the spherical model with the density contrast 
calculated by the linear theory. We define the critical density contrast 5 C that is the value, extrapo- 
lated to the present time (t = to) using linear theory, of the overdensity which collapses at t = t co \\ 
in the exact spherical model. It is given by 



<5 c (icoii) = \D{t Q ) 



V^coll/ 



(fio < 1) (5) 



(Lacey and Cole 1993), where D(t) is the linear growth factor given by equation (A13) of Lacey 
and Cole (1993) and t n = 7 riJ^ 1 fi (l - ^o)^ 3/2 - 

For a power-law initial fluctuation spectrum P(k) ex k n , the rms amplitude of the linear mass 
fluctuations in a sphere containing an average mass M at a given time is 5 a M-(™+ 3 )/ 6 . Thus, 
the virial mass of clusters which collapse at t co n is related to that at to as 



M vir (i co u) = M virj0 



6 c (t 



collj 



L <*c(to) 



-6/(n+3) 

(7) 



Here, M v i r) o(= M v i r [to]) is regarded as a variable because actual amplitude of initial fluctuations 
has a distribution. We relate t = t co u to the collapse or formation redshift z co \\, which depends on 
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cosmological parameters. Thus, M v i r is a function of z co \\ as well as M v ; ri o. This means that for 
a given mass scale M v ; r , the amplitude of initial fluctuations takes a range of value, and spheres 
containing a mass of M v ; r collapse at a range of redshift. In the following, the slope of the spectrum 
is fixed at n = — 1. It is typical of the scenario of standard cold dark matter for a cluster mass 
range, and is consistent with observations as shown in Paper II. 



The virial radius and temperature of a cluster are then calculated by 

'3M V - xV3 



4irp v 



(8) 



rp _ ^ra H GM vir ^ 

3/CB r vir 

where fj,(= 0.6) is the mean molecular weight, mu is the hydrogen mass, k-g is the Boltzmann 
constant, G is the gravitational constant, and 7 is a fudge factor which typically ranges between 1 
and 1.5. In Paper II, we adopted the value 7 = 1. Note that the value of 7 is applied only to dark 
matter, but not to gas, because we do not assume Tg as — T v ; r here. We emphasize that M v [ r , p v ; r , 
and r v i r are the virial mass, density, and radius at the time of the cluster collapse, respectively. 



2.2. Shocks and Hydrostatic Equilibrium 

To study the effect of preheating, we here adopt a very simple model as a first step. When 
a cluster collapses, we expect that a shock wave forms and the infalling gas is heated. In order 
to derive the postshock temperature, we use a shock model of Cavaliere et al. (1998). For a 
given preshock temperature Ti, the postshock temperature T2 can be calculated from the Rankine- 
Hugoniot relations. Assuming that the shock is strong and that the shock front forms in the vicinity 
of r vir , it is approximately given by 

kB T 2 = -fel + h B T, (10) 

(Cavaliere et al. 1998), where (f>(r) is the potential at r. According to the virial theorem and the 
continuity when Ti approaches zero, we should take — </>(r v i r )/3 = k^T vu . For r < r v i r , we assume 
that the gas is isothermal and hydrostatic, and that the matter accretion after the cluster collapse 
does not much change the structure of the central region of the cluster significantly, as confirmed 
by numerical simulations (e.g. Takizawa and Mineshige 1998). It is to be noted that even if the 
density profile of dark matter is represented by the universal profile (Navarro et al. 1996, 1997), 
it is not inconsistent with the isothermal (3 model of gas represented by equation (2) (Makino et 
al. 1998; Eke et al. 1998) within the present observational scopes. On these assumptions, the gas 
temperature in the inner region of a cluster is T gas = T2 , and the mass within r of the cluster center 
Mdm is related to the density profile of intracluster gas, p gas , by 

M = ~~pm^G r ~~dlirr~ < U > 
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Since MDM(»"vir) = M v - ir , equations (9) and (11) yield 

(12) 



rp __1 T dlnp ga 

-*- vir — ~ -*- c; 



3 ° dlnr 

Defining /? = T v i r /T gas , the gas density profile is thus given by 

PgasW oc r-^h , (13) 

as long as (dlnp gas /dln r) is nearly constant. 

Equation (10) shows that in this model (3 is a function of only T v ; r when Ti is regarded as an 
external parameter, that is, 

p = T ; ir ; - . (14) 

P T vir + (3/2)Ti 1 ' 

Since T v ; r = T gas f3, equation (14) is written as 

g= T 8M (3/m (15) 

-tgas 

Thus, the /3 — T gas relation can be used to determine T\ by comparing with the observation. Since 
both T v i r and r v ; r are the two-parameter families of z co \\ and M v i r) o (equations [3], [7], [8] and [9]), 
equation (14) shows that (3 can be represented by r v ; r as (3 = /3(r v [ T , M v ; r) o), if T\ is specified. 
Recent numerical simulations suggest that the structure of central region of clusters is related 
to z co n (Navarro et al. 1997), and in particular tdm,c is proportional to r v - lv (Salvador-Sole et al. 
1998; Makino et al. 1998). Therefore, if we assume that rrjM,c = r c an d that r vir /r c is constant 
as in Paper II, T\ can also be determined by comparing the theoretical prediction of the (3 — r c 
relation with the observation. Since a spherical collapse model predicts r v i r (z co n = 0) ~ 4 Mpc and 
observations show that r c (z co n = 0) ~ 0.5 Mpc (Figure lb in Paper II), we adopt r v i r /r c = 8 from 
now on. Thus, we obtain (3 = 0(8r c [z co \\, M^o], M v i r; o). 



3. Results and Discussion 

3.1. (3 — T gas and [3 — r c Relations 

Using the model constructed in §2.2, we predict relations between (3 and T gas , and between (3 
and r c . 

If T\ is mainly determined by the energetic winds generated in the forming galaxies or quasars 
before the formation of clusters, T\ should be constant if subsequent adiabatic heating or cooling 
is neglected. However, if, besides the winds, the gravitational energy released in the subclusters, 
which later merged into the observed clusters, contributes to Ti, we expect that Ti has a distribution 
produced by different merging histories. In order to determine the distribution in detail, we must 
calculate the merging histories by Monte Carlo realizations as Cavaliere et al. (1997, 1998) did. In 
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this study, however, we consider the scatter by investigating a range of T\ for simplicity We show 
in Figure 1 the (3 — T gas relation for T\ = 0.5, 1, and 2 keV. The observational data are overlaid. 
Since equation (2) is approximated to be /9 gas (?*) oc r~ 3/3obs for r >> r c , the relation 

(3 = 7/3obs (16) 

is obtained by comparing the relation (13). Thus, in the following figures, the observed values of 
Pobs are converted by relation (16). In Figure 1 we assumed 7 = 1. As the data, we use only 
relatively hot (T gas > 3 keV) and low redshift (z < 0.1) clusters obtained by Mohr et al. (1999) and 
Peres et al. (1998). Instead of f3 hs, Peres et al. (1998) present velocity dispersions corresponding 
to gravitational potential well, 0"dcproj 5 derived with the deprojection method, ignoring velocity 
dispersion anisotropies and gradients . Thus, for the data we assume that k^T^u = ^"iH0"d epro j and 
define (3 as T vir /T gas . Figure 1 shows that the observational data are consistent with 0.5 < Ti < 2 
keV but it seems that a single value of T\ does not represent the range of data. The preheating 
(Ti > 0) is expected to reduce (3 of the clusters with small T gas (Figure 1). At first glance, no 
correlations between (3 and T gas are recognized observationally in this temperature range. However, 
some reports on the existence of a weak correlation have been made when clusters with lower T gas are 
included (e.g. Horner et al. 1999). Thus, our prediction is not inconsistent with the observations. 

As discussed in §2.2, the (3 — r c relation is represented by two parameters z co \\ and M v - ir ^, for 
a given value of T\. The results are shown in Figure 2 for 7 = 1. Figure 2a and 2b are for = 1 
and 0.2, respectively. For comparison, we also present observational data (Mohr et al. 1999; Peres 
et al. 1998). As was in Paper I, for the data of Mohr et al. (1999) we use here only the component 
of surface brightness reflecting the global structure of clusters, although the central component 
(so-called cooling flow component) may also have formed in the scenario of hierarchical clustering 
(Fujita and Takahara 2000). 

The mass M v i rj o corresponds to the mass of clusters collapsed at z ~ and takes a range of value 
due to the dispersion of initial density fluctuation of the universe. Since observations and numerical 
simulations show M v i rj o ~ 10 15 M Q (Evrard et al. 1996), the observational data are expected to lie 
between the two lines of M virfi = 5 x 10 14 M (arc BC) and M virfi = 5 x 10 15 M (arc AD) for fixed 
T\ in Figure 2. Note that the distribution of M v i rj o degenerates on the lines in Figure 1. In Figure 2, 
the positions along the arcs AD and BC indicate the formation redshifts of the clusters. When 
fio = 1, most of the observed clusters should have collapsed at z ~ because clusters continue 
growing even at z = (Peebles 1980). Thus, the cluster data are expected to be distributed along 
the part of the lines close to the point of z co \\ = (segment AB). In fact, calculations done by 
Lacey and Cole (1993), and Kitayama and Suto (1996) show that if S^o = 1, most of present day 
clusters (M v j r ~ 1O 14_15 M ) should have formed in the range of z co \\ < 0.5 (parallelogram ABCD in 
Figure 2a). In contrast, when Qq = 0.2, the growth rate of clusters decreases and cluster formation 
gradually ceases at z < l/^o ~~ 1 (Peebles 1980). Thus, in Figure 2b, cluster data are expected to 
be distributed between the points of z co n = (segment AB) and z co n = — 1 (segment CD) 

and should have a two-dimensional distribution (parallelogram ABCD). Thus, compared with the 
observations, the models in Figure 2 show that T\ ~ 1 keV and Oq < 1 are preferred. The latter 
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result is quite consistent with that of Paper II, where we found that the T gas — r c relation suggests 
0,q < 1. Since (3 is related to T gas by equation (15), the [3 — r c relation is equivalent to the T gas — r c 
relation for a fixed value of T±. Note that in Figure 2 predicted regions corresponding to different 
T\ overlap each other; this implies that the position of a cluster in Figure 2 does not uniquely 
correspond to Ti in contrast to Figure 1. For a given (3 and r c , larger T\ corresponds to larger 
-^vir,o or larger amplitude of the initial fluctuation. 

The dispersion of Ti appears to be caused by gravitational heating in subclusters that are to 
merge to the cluster (T gas > 3 keV) at the time of the cluster formation. In fact, Figure 2b shows 
that observed clusters are situated close to the line of M v ; ri o = 5 x 10 15 M (arc AD) when T\ ~ 2 
keV, while they are situated close to the line of M N - 1T ^ = 5 x 10 14 M (arc BC) when 0.5 < T\ < 1 
keV. Moreover, Figure 1 suggests that clusters with large T gas favor large T\. These may reflect 
that clusters with larger (smaller) M v i rj o or T gas tend to have more (less) massive progenitors with 
larger (smaller) Ti, although these are only loose tendencies, and we need more samples and more 
improved models to obtain a definite conclusion. 

Note that gravitational heating in subclusters itself is a self-similar process and does not 
modify self-similar scaling relations such as the luminosity-temperature relation (e.g. Eke et al. 
1998). Thus, an additional entropy other than expected from purely gravitational assembly of a 
cluster must be injected into the gas. Valageas and Silk (1999) investigate the entropy evolution of 
intergalactic medium (IGM) and find that clusters with T v i r ~ 0.5 keV are affected by the additional 
entropy when it is generated by quasar heating. This is because the additional entropy is comparable 
to the entropy generated by gravitational collapse of the clusters. In other words, the adiabatic 
compression of the gas from the preheated IGM alone can heat the gas up to T a d jC i ~ 0.5 keV. 
Therefore, in addition to the gravitational processes in subclusters, the preheating may significantly 
contribute to Ti, and the lower bound of which is given by T a( j c i. If T a d )C i ~ 0.5 keV, this is consistent 
with our result (Figures 1 and 2). Valageas and Silk (1999) also investigate the case when only 
supernova heating is taken into account and quasar heating is ignored. The result is T^ci < 0.1 
keV. In this case, and effects of the preheating is small and we expect that (3 does not much depend 
on T gas and r c , although (3 would have a scatter owing to the difference of merging history. This is 
inconsistent with the observations. The insufficient power of the supernova heating is also suggested 
by Wu et al. (1999) (but see Loewenstein 1999). Another possible source of heating is that due 
to shocks forming at higher redshift on the largest scales, such as filaments and sheets. Cen and 
Ostriker (1999) indicate that most of baryons at low redshift should have a temperature in the 
range of 10 5 — 10 7 K. The relatively large value of T\ may reflect this temperature. 

We also investigate the case of 7 = 1.2 and VLq = 0.2, which are presented in Figure 3. In this 
case, the model of Ti = 0.5 keV is preferred especially for the data obtained by Mohr et al. (1999). 
This means that 7 and Ti are correlated and they cannot be determined independently. However, 
the model of 7 > 1.2 is inappropriate because (3 = 7/3 bs exceeds unity for some observational 
data while relation (14) or (15) limits [3 to less than one. If a cluster is not isothermal, the 
temperature in the central region T gas should be larger than T2 (Cavaliere et al. 1998). In this case, 
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the discrepancy between the model and the observations is more significant. Thus, it seems to be 
difficult to construct a model that predicts (3 > 1. 



3.2. The Fundamental Band and Plane 

It is interesting to investigate whether the gas distribution in clusters derived above is consistent 
with the observations of central gas fraction, and the fundamental band and plane we found in 
Paper I. The shapes of the band and plane are also related to the origin of the observed relation of 
Lx oc Tg as (Paper I). We did not explore the origin of the variation of the central gas mass fraction 
in previous papers, where /3 G bs was regarded as constant. Below, we will show that this is related 
to the variation of (3. 

From relation (13), the gas density at the cluster core is approximately given by 

Pgasfi = /Vs( r vir)(7vir/V c ) 3/3//7 , (17) 

where r v j r and P(T v i r , T±) are functions of z co \\ and M v i rj o (§2), and p g as{r) is the gas density at radius 
r from the cluster center. We assume that the profile of dark matter is isothermal (pom °c r ~ 2 ) 
at least for r c < r < r v - 1T , and pr>M,c = 64/9 v ; r . Moreover, we assume that the average gas fraction 
within radius r v - 1T is f gas {r v i T ) = 0.25 regardless of z co ii and M v ; rj o. The value of / gas is nearly the 
largest gas mass fraction of observed clusters (e.g. David et al. 1995; Ettori and Fabian 1999). On 
these assumptions, the central gas density and the gas fraction at the cluster core are respectively 
given by 

( P\ fr vir \ 3/3/7 

Pgas,0 = ( 1 - — J /gas Pvir(^coll) ( ~ ) 

= 0.25 (l - ^ Pvi r(z coll ) 8 3 ^ (18) 

and 

/ g as(0) = 0.25 (l - ?\ 8*W^ , (19) 

where / ga s(0) = /O g as,o//ODM,c is the gas fraction at the cluster center. The above equations are valid 
when (3 < 7. Note that in Paper II, we derive the central gas density according to the relation 
Pgas,o oc Pvir/gas(0), in which / gas (0) is differently determined by the observations 3 . In contrast, in 
equation (18), we derive p gas ,o assuming that / gas (^vir) = constant. 

The above model values (equation [18] or [19]) can be obtained from observational data. Using 
equations (8), (9), and (16) we obtain 



ATrGfimu (8r c ) 2 ' 



(20) 



3 In Papers I and II, we assumed that T gas = T v ; r and did not take account of the variation of f3 b B when we derive 
/gas(0) from observations of p gas ,o, r c , and T gas . 
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where we used the relations of T v ; r = f3T gas and r c = r v ; r /8. Thus, using equation (16), the right 
hand of equation (18) should be written as 



model _ n O r/o /-i o \ ^jj-g^gas 8 3/3 ° bs 

p gaSi0 = 0.25/3 obs (l - /Jobs)^^ (8rc)2 • (21) 



Hence, p™* can be derived from the observable quantities r c , T gas , and /? Q bs- Figure 4 displays a 
plot of p g as,o and P^a° d o el based on the data obtained by Mohr et al. (1999). Note that Peres et al. 
(1998) do not present /9 ga s,o- We do not show the uncertainties of p^ s d Q l to avoid complexity. Here 
we use only /9 gas ,o corresponding to the global cluster component as we did in Paper I. 

Figure 4 shows that p ga s,o we h agrees with Pg^o' although /5g as ,o is slightly smaller than /9^ a ° d o 1 
for clusters with large p™^^ 1 - Thus, we conclude that the variation of / gas (0) is due to that of the 
slope parameter of the gas distribution (3 within r v ; r . One possible reason for the slight disagreement 
between /5™° d el and p gas ,o is an uncertainty of the value of ig as (?Vir)- Another is the influence of 
central excess emission of clusters. When the distance to a cluster is relatively large, the center and 
global surface brightness components may not be distinguished even if the two components exist. 
In this case, the cluster may be considered that it has only a global component. However, when 
the central emission is strong, the fitting of the surface brightness profile by one component may 
be affected by the central emission and may give a smaller core radius than the real. This may 
make Pg^° d o el large for the clusters. In fact, clusters with /0g^° d el > 3 x 1CT 26 g cm -3 are regarded 
by Mohr et al. (1999) as having only one (global) component of surface brightness. Note that core 
radii derived by Peres et al. (1998) may be less affected by the central emission because they take 
account of cooling flows and the gravitation of central cluster galaxies, which are responsible for 
the central emission, for all clusters they investigate (Figures 2 and 5b). 

We present the theoretically predicted relations among p gas ,0) r vin and T gas in Figure 5. 
Although these relations are presented in Paper II using the observed relation between / gas (0) 
and Mdm, c , here we plot them by directly using (3. For lines in Figure 5, we use the rela- 
tion T gas = T v i T /(3. For comparison, we plot the observational data in the catalogue of Mohr 
et al. (1999) and Peres et al. (1998). For the data, we use the relation r v - 1T = 8r c . Figure 5 
shows that our model can well reproduce the band distribution of observational data in the 
(/0 g as,O) r o T gas )-space. Moreover, our model can explain the planar distribution of the observa- 
tional data. In Paper I, we find that the observational data satisfy the relation of the fundamental 
plane, p g alVc' 65T gas' 6 ° a constant. For M Vufi ~ 10 15 M and z coU < 2, our model with S7 = 0.2 
and T\ = 1 keV predicts the plane of /fg a s^o r c' 64 ^gas' 7 ° ^ constant, which is approximately con- 
sistent with the observation. Note that the index of /0 ga s,o is somewhat smaller than the observed 
value considering the uncertainty (~ 0.1), which may be related to the slight disagreement between 
Pgaslj 1 and p g as,o (Figure 4). The plane is represented by the two parameters, z co \\ and M, 



vir,0) 



as discussed in §2. Since the cross section of the fundamental plane corresponds to the observed 
Lx~ Tgas relation, and the fundamental plane corresponds to the observed dependence of / gas (0) on 
Pdm,c and Mdm,c (Paper I), our model can also reproduce these relations. These results strengthen 
our interpretation that the difference of gas distribution among clusters is caused by heating of the 



- 11 - 



gas before the cluster collapse and by shock heating at the time of the cluster collapse (equation 
[14] or [15]). 

4. Conclusions 

We have investigated the influence of heating before cluster collapse and shocks during cluster 
formation on the gas distribution in the central region of clusters of galaxies. We assumed that the 
core structure has not much changed since the formation of a cluster. Using a spherical collapse 
model of a dark halo and a simple shock model, we predict the relations among the slope of gas 
distribution (3, the gas temperature T gas , and the core radius r c of clusters. By comparing them with 
observations of relatively hot (> 3 keV) and low redshift clusters, we find that the temperature 
of the preheated gas collapsed into the clusters is about 0.5 — 2 keV. Since the temperature is 
higher than that predicted by a preheating model of supernovae, it may reflect the heating by 
quasars or gravitational heating on the largest scales at high redshift. Moreover, gravitational 
heating in subclusters assembled when the clusters formed also seems to affect the temperature 
of the preheated gas and produce the dispersion in the preheating temperature. Assuming that 
the global gas mass fraction of clusters are constant, we predict that the gas mass fraction in the 
core region of clusters should vary correlating with (5 through a simple law, which is shown to be 
consistent with the observations. Thus, we conclude that the variation of the gas mass fraction in 
the cluster core is due to the shock heating of preheated gas. Furthermore, we have confirmed that 
the observed fundamental plane and band of clusters are reproduced by the model even when the 
effects of preheating are taken into account. Thus, major conclusions about the cluster formation 
and cosmology obtained in our previous papers are not changed. 

We thank for A. C. Edge, C. S. Frenk, and T. Kodama for useful discussions. This work was 
supported in part by the JSPS Research Fellowship for Young Scientists. 
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Fig. 1. — The relation between (5 and T gas when 7 = 1. Dashed line: T\ = 0.5 keV. Solid line: 
Ti = 1 keV. Dotted line: T x = 2 keV. The observational data obtained by Mohr et al. (1999) (filled 
circles) and Peres et al. (1998) (open squares) are overlaid. 
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Fig. 2. — Theoretical predictions of (3 — r c relation in the case of (a) f2o = 1-0 and 7=1 (b) 
tt = 0.2 and 7 = 1. Dashed line: T x = 0.5 keV. Solid line: T x = 1 keV. Dotted line: Ti = 2 keV. 
The arcs AD and BC correspond to M virj0 = 5 x 10 15 M and 5 x 10 14 M , respectively. The 
segments AB corresponds to z co \\ = and the segments CD corresponds to z co n = 0.5 in (a) and 
z co ii = 4 in (b). The observational data obtained by Mohr et al. (1999) (filled circles) and Peres et 
al. (1998) (open squares) are overlaid. 
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Fig. 3. — The same as Figure 2b but for 7 = 1.2 




Fig. 4. — The relation between /J^° d cl of equation (21) and p gas fi. The line corresponds to /9 gas ,o = 

„modcl 
^gas,0 
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Fig. 5. — Theoretical predictions of (a) radius-density relation (b) radius-temperature relation (c) 
density-temperature relation when 7 = 1. Solid line: f^o = 0.2 and M v i ri o = 5 x 10 15 M . Dotted 
line: O = 0.2 and M Virfi = 5 x 10 14 M . Dashed line: n = 1.0 and M viri0 = 5 x 10 15 M . 
Dash-dotted line: O = 1-0 and M virfi = 5 x 10 14 M . The open diamonds, triangles, and circles 
correspond to the collapse redshifts of z co \\ = 0, z co \\ = 0.5, and z co \\ = 4, respectively. We assume 
that r v ; r = 8r c . The observational data obtained by Mohr et al. (1999) (filled circles) and Peres et 
al. (1998) (open squares, only in Figure 5b) are overlaid. 
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